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Abstract: We review the global issues associated to gauge fixing ambiguities 
and their consequence for glueball spectroscopy. To avoid infrared singulari- 
O I ties the theory is formulated in a finite volume. The examples of a cubic and 

^ ■ spherical geometry will be discussed in some detail. Our methods are not 

' powerful enough to study the infinite volume limit, but the results clearly in- 

dicate that for low-lying states, wave functionals are sensitive to global gauge 
copies which we will argue is equivalent to saying that they are sensitive to 
^ I the geometric and topological features of configuration space. 

1 Introduction 

in . 

G^ ■ Gauge theories with non-Abelian gauge groups lead not only to more complicated poten- 
^ [ tials, but also to complicated kinetic terms in the Lagrangian or Hamiltonian. The latter 
are a manifestation of the non-trivial Riemannian geometry of the physical configuration 
space 0, formed by the set of gauge orbits A/Q {A is the collection of connections, Q 
^ ■ the group of local gauge transformations). Most frequently, coordinates of this orbit space 
are chosen by picking a representative gauge field on the orbit in a smooth and prefer- 
^ [ ably unique way. It is by now well known that linear gauge conditions like the Landau or 
^ I Coulomb gauge suffer from Gribov ambiguities 0. The reason behind this is that topo- 
logical obstructions prevent one from introducing afiine coordinates in a global way. In 
principle therefore, one can introduce different coordinate patches with transition functions 
to circumvent this problem |Q. One way to make this specific is to base the coordinate 
patches on the choice of a background gauge condition. One could envisage to associate 
to each coordinate patch ghost fields and extend ERST symmetry to include fields with 
non-trivial "Grassmannian sections", although such a formulation is still in its infancy. 
Interesting conjectures concerning non-perturbative spontaneous breakdown of BRST in- 
variance are implied by the work of Zwanziger and collaborators ^, but will not be 
discussed here. 

We will pursue, however, the issue of finding a fundamental domain for non-Abelian 
gauge theories and its consequence for the glueball spectrum in intermediate volumes. 



^Talk at the ECT* workshop "Non-perturbative approaches to QCD", July 10 - 29, 1995, Trento, Italy 
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The finite volume context allows us to make reliable statements on the non-perturbative 
contributions, because asymptotic freedom guarantees that at small volumes the effective 
coupling constant is small, such that high-momentum states can be treated perturba- 
tively. Only the lowest (typically zero or near-zero momentum) states will be affected 
by non-perturbative corrections. We emphasize that it is essential that gauge invariance 
is implemented properly at all stages. We will describe the results mainly in the con- 
text of a Hamiltonian picture with wave functionals on configuration space. Although 
rather cumbersome from a perturbative point of view, where the covariant path integral 
approach of Feynman is vastly superior, it provides more intuition on how to deal with 
non-perturbative contributions to observables that do not vanish in perturbation theory. 
An essential feature of the non-perturbative behaviour is that the wave functional spreads 
out in configuration space to become sensitive to its non-trivial geometry. If wave function- 
als are localized within regions much smaller than the inverse curvature of the field space, 
the curvature has no effect on the wave functionals. At the other extreme, if the config- 
uration space has non-contractible circles, the wave functionals are drastically affected by 
the geometry, or topology, when their support extends over the entire circle. Instantons 
are of course the most important examples of this. Not only the vacuum energy is affected 
by these instantons, but also the low-lying glueball states and this is what we are after to 
describe accurately, albeit in sufficiently small volumes. The geometry of the finite volume, 
to be considered here, is the one of a three-torus |rU| and a three-sphere [|TT], 0. 

Some complications with so-called reducible connections (connections with a non- 
trivial stabilizer, i.e. subgroup of the gauge group that leaves the connection invariant). 



discussed in detail in ref. |TJ], will not occur in the sector where we study the dynamics. 
Although A = is such a reducible connection, and would give rise to a curvature singu- 
larity in configuration space, we know perfectly well how to deal with it by not fixing the 
constant gauge transformations, which form the stabilizer of this vector potential. Indeed 
the Coulomb gauge does not fix this gauge degree of freedom. We simply demand that the 
wave functional is in the singlet representation under the constant gauge transformations. 

Finally we should mention that recently the issue of summing over all gauge copies ||T^ 
has been revived by studying a simple soluble "gauge" model. It is shown in this 
model that singularities which arise through the vanishing of the Faddeev-Popov deter- 
minant, which is associated to the Gribov horizon, can be properly accounted for. It is 
not clear that this can also be achieved in the case of Yang-Mills theories, where things 
are far less tractable, although a recent paper by Fujikawa gives arguments in favour 
of this on the basis of a careful BRST analysis. It is argued in ref. [|16| that selecting 
only one gauge copy requires imposing artificial boundary conditions, which changes the 
topology of the problem. This is exactly contrary to what was this authors original mo- 
tivation, namely that the non-trivial topology enforces one to consider those boundary 
conditions |T^. In principle, as we will see, these boundary conditions are well defined 



and only involve identifications by a gauge transformation, under which we know exactly 
how the wave functional transforms. Although it is correct to state that it is difficult to 
develop the Feynman rules for perturbative calculations, it should be emphasized that the 
global issues that arise from the gauge invariance of the theory, leading to a topologically 
and geometrically non-trivial configuration space, give rise to non-perturbative effects that 
are not expected to be predicted reliably by perturbation theory. The issue is to isolate 
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the relevant non-perturbative contributions and to include them in a way that will not 
violate the gauge invariance of the theory. In that sense one can even take the extreme 
point of view in ref. [jl9[ where it is attempted to eliminate all gauge degrees of freedom 



by formulating the theory in terms of electric or magnetic fields, although this leads to all 
sorts of technical difficulties, but it has the advantage that approximations made after this 
reformulation is implemented do not break the gauge invariance. 

In the following we provide a quick review of the definition of the fundamental domain. 
For gauge theories on a torus, when restricting oneself to the zero-momentum modes, 
the fundamental domain is simple and leads to a reliable way of computing the low-lying 
glueball spectrum in volumes up to half a cubic fermi. We then describe the situation for 
the three sphere, where also a rather complete picture of the fundamental domain for the 
low-lying modes has been achieved. 

2 Gribov and fundamental regions 

An (almost) unique representative of the gauge orbit is found by minimizing the norm 
of the vector potential along the gauge orbit [0, 

FA{g) = rAf = - J^^d'x tr l^(g-'Ag + g-%g 

where the vector potential is taken anti-hermitian. For SU{2), in terms of the Pauli 
matrices r^, one has: 

g{x) = exp{X{x)), X{x)=zX'^{x)^. (2) 
Expanding around the minimum of eq. (P, one easily finds: 
= + ti {Xd^Ai) + f ti {X^FP{A)X) 

JM JM 

+\f tr (Xp„X],9.X]) + i- / tr ([AX,X][9.X,X]) + 0(X^). (3) 

O JM iZ JM ^ ' 

Where FP(y4) is the Faddeev-Popov operator (ad(yl)X = [y4,X]) 

FP{A) = -dMA) = -d^-dM{Ai). (4) 

At any local minimum the vector potential is therefore transverse, diAi = 0, and 
FP{A) is a positive operator. The set of all these vector potentials is by definition the 
Gribov region Q. Using the fact that FP{A) is linear in A, Q is seen to be a convex subspace 
of the set of transverse connections P. Its boundary dQ is called the Gribov horizon. At 
the Gribov horizon, the lowest eigenvalue of the Faddeev-Popov operator vanishes, and 
points on dQ are hence associated with coordinate singularities. Any point on dQ can be 
seen to have a finite distance to the origin of field space and in some cases even uniform 
bounds can be derived ETl, . 
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The Gribov region is the set of local minima of the norm functional (3) and needs to 
be further restricted to the absolute minima to form a fundamental domain, which will be 
denoted by A. The fundamental domain is clearly contained within the Gribov region. To 
show that also A is convex we note that 

VAr-Ur = / tr (A^) -/ tr [{g^'Ag + g-'d,g)') 

= JtT (g^FP._ {A) g) ^ {g, FP. {A) g) , (5) 

where FPi{A) is the Faddeev- Popov operator generalized to the fundamental representa- 
tion. Or for the gauge group SU(2) we have 

FP,{A) = -dl - \A1T^d.. (6) 

Here Tf are the hermitian gauge generators in the spin-t representation: 

T| = |, = ad(|). (7) 

They are angular momentum operators that satisfy = t{t + l)id. At the critical points 
y4 G r of the norm functional, (recall F = {A G = 0}), FPt{A) is an hermitian 

operator. Furthermore, FPi{A) in that case coincides with the Faddeev- Popov operator 
FP{A) in eq. (P. We can define A in terms of the absolute minima (apart from the 
boundary identifications) over 5^ G ^ of (^g, FPi{A) g^ 

A = {A G F| min (g, FP. {A) g) = 0}. (8) 

Using that FPi{A) is linear in A and assuming that A^^^ and A^'^'^ are in A and therefore 

satisfy the equation min^^^ (^g,FPi{A) g^ = 0, we find that A = sA^^'> + (1 — s)A^'^'^ 
satisfies the same identity for all s G [0, 1] (such that both s and (1 — s) are positive). The 
line connecting two points in A, therefore lies within A. 

Its interior is devoid of gauge copies, whereas its boundary dA will in general contain 
gauge copies, which are associated to those vector potentials where the absolute minima of 
the norm functional are degenerate [0 . If this degeneracy is continuous one necessarily has 
at least one zero eigenvalue for FP{A) and the Gribov horizon will touch the boundary of 
the fundamental domain at these so-called singular boundary points. By singular we mean 
here a coordinate singularity. We sketch the situation in figure 1. As mentioned already in 
the introduction, the constant gauge degree of freedom is not fixed by the Coulomb gauge 
condition and therefore one still needs to divide by G to get the proper identification 

A/G = A/g. (9) 

Here A is considered to be the set of absolute minima modulo the boundary identifica- 
tions, that remove the degenerate absolute minimum. It is these boundary identifications 
that restore the non-trivial topology of A/ Q. Furthermore, the existence of non-contractible 
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spheres [§| allows one to argue that singular boundary points are to be expected |T8|. Con- 



sider the intersection of a d-dimensional non-contractible sphere with A. The part in the 
interior of A is contractible and it can only become non-contractible through the bound- 
ary identifications. It was falsely stated in ref. |jl8| that this can only be the case if for 
the ((i-l)-dimensional intersection with the boundary, all points are to be identified. As 
this would imply degeneracy of the the norm functional on a ((i-l)-dimensional subspace, 
there would be at least d-1 zero-modes for the Faddeev-Popov operator. Although each 
connected component of the intersection with the interior of A is contractible, there can be 
more than one connected component. In the case of two such components one can make a 
non-contractible sphere by identifing points of the ((i-l)-dimensional boundary intersection 
of the first connected component with that of the second and there is no necessity for a 



continuous degeneracy 23 



Not all singular boundary points, even those associated with continuous degeneracies, 
need to be associated with non-contractible spheres. Note that absolute minima of the 
norm functional are degenerate along the constant gauge transformations, this is a trivial 
degeneracy, also giving rise to trivial zero-modes for the Faddeev-Popov operator, which 
we ignore. The action of G is essential to remove the curvature singularities mentioned 
above and also greatly facilitates the standard Hamiltonian formulation of the theory 
There is no problem in dividing out G by demanding wave functionals to be gauge singlets 
(colourless states) with respect to G. In practice this means effectively that one minimizes 
the norm functional over Q/G. When a singular boundary point is not associated to a 
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continuous degeneracy, the norm functional undergoes a bifurcation, when we move from 
inside to outside the fundamental (and Gribov) region. The stable absolute minimum 
turns into a saddle point and two stable minima appear, as indicated in figure 2. These are 
necessarily gauge copies of each other. The gauge transformation is homotopically trivial 
as it reduces to the identity at the bifurcation point, evolving continuously from there on. 
The explicit example of ref. is one where the connection is reducible and it can be 
shown that the two stable minima that appear after the bifurcation are gauge copies by 
constant gauge transformations |24], the situation being subtle as the gauge transformation 
9i^92 (see figure 2) is in this case not constant, but it is in the stabilizer of the appropriate 
gauge field, up to a constant gauge transformation. Examples of bifurcations at irreducible 
connections were explicitly found for 5*^, see ref. (app. A). We will come back to this. 

Also Gribov's original arguments for the existence of gauge copies (showing that 
points just outside the horizon are gauge copies of points just inside) can be easily under- 
stood from the perspective of bifurcations in the norm functional. It describes the generic 
case where the zero-mode of the Faddeev-Popov operator arises because of the coalescence 
of a stable minimum with a saddle point with only one unstable direction, see ref. [0 for 
more details and a discussion of the Morse theory aspects that simplify the bifurcation 
analysis. 
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3 Gauge fields on the three-torus 



Homotopical non-trivial gauge transformations are in one to one correspondence with non- 
contractible loops in configuration space, which give rise to conserved quantum numbers. 
The quantum numbers are like the Bloch momenta in a periodic potential and label rep- 
resentations of the homotopy group of gauge transformations. On the fundamental do- 
main the non-contractible loops arise through identifications of boundary points (as will 
be demonstrated quite explicitly for the torus in the zero- momentum sector). Although 
slightly more hidden, the fundamental domain will therefore contain all the information 
relevant for the topological quantum numbers. Sufficiently accurate knowledge of the 
boundary identifications will allow for an efficient and natural projection on the various 
superselection sectors (i.e. by choosing the appropriate "Bloch wave functionals"). All 
these features were at the heart of the finite volume analysis on the torus [|1^] and we 
see that they can in principle be naturally extended to the full theory, thereby including 
the desired 6 dependence. In the next section this will be discussed in the context of 
the three-sphere. In ref. we proposed formulating the Hamiltonian theory on coordi- 
nate patches, with homotopically non-trivial gauge transformations as transition functions. 
Working with boundary conditions on the boundary of the fundamental domain is easily 
seen to be equivalent and conceptually much simpler to formulate. If there would be no 
singular boundary points this would have provided a Hamiltonian formulation where all 
topologically non-trivial information can be encoded in the boundary conditions. Still, 
for the low-lying states in a finite volume, both on the three-torus and the three-sphere, 
singular boundary points will not play an important role in intermediate volumes. 

Probably the most simple example to illustrate the relevance of the fundamental do- 
main is provided by gauge fields on the torus in the abelian zero-momentum sector. For 
definiteness let us take G = SU (2) and Ai = i^iT^ {L is the size of the torus). These modes 
are dynamically motivated as they form the set of gauge fields on which the classical po- 
tential vanishes. It is called the vacuum valley (sometimes also referred to as toron valley) 
and one can attempt to perform a Born-Oppenheimer-like approximation for deriving an 
effective Hamiltonian in terms of these "slow" degrees of freedom. To find the Gribov 
horizon, one easily verifies that the part of the spectrum for FP{A) that depends on C, is 
given by Xf{C) = 27rn • (27rn ± C), with n ^ an integer vector. The lowest eigenvalue 
therefore vanishes if = ±27r. The Gribov region is therefore a cube with sides of length 
in, centred at the origin, specified by \Ck\ < 2n for all k, see figure 3. 

The gauge transformation g(^k) = ewi'^^'i^kTs/ L) maps Ck to Ck + 2n, leaving the other 
components of C untouched. As g(^k) is anti-periodic it is homotopically non-trivial (they 
are 't Hooft's twisted gauge transformations [^). We thus see explicitly that gauge copies 
occur inside fl, but furthermore the naive vacuum A = has (many) gauge copies under 
these shifts of 27r that lie on the Gribov horizon. It can actually be shown for the Coulomb 
gauge that for any three-manifold, any Gribov copy by a homotopically non-trivial gauge 
transformation of A = will have vanishing Faddeev-Popov determinant |18|. Taking the 
symmetry under homotopically non-trivial gauge transformations properly into account is 
crucial for describing the non-perturbative dynamics and one sees that the singularity of 
the Hamiltonian at Gribov copies of A = 0, where the wave functionals are in a sense 
maximal, could form a severe obstacle in obtaining reliable results. 
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Figure 3: A two dimensional slice of the vacuum valley along the (Ci,C2) plane. The fat 
square give the Gribov horizon, the grey square is the fundamental domain. The dots at 
the Gribov horizon are Gribov copies of the origin. 



To find the boundary of the fundamental domain we note that the gauge copies C = 
(vr, C2, C3) and C = (— vr, C2, C3) have equal norm. The boundary of the fundamental do- 
main, restricted to the vacuum valley formed by the abelian zero-momentum gauge fields, 
therefore occurs where \Ck\ = vr, well inside the Gribov region, see figure 3. The bound- 
ary identifications are by the homotopically non-trivial gauge transformations g(^k)- The 
fundamental domain, described by \Ck\ < vr, with all boundary points regular, has the 
topology of a torus. To be more precise, as the remnant of the constant gauge transforma- 
tions (the Weyl group) changes C to —C, the fundamental domain A/G restricted to the 
abelian constant modes is the orbifold T^/Z2. Generalizations to arbitrary gauge groups 
were considered in ref. (The fundamental domain turns out to coincide with the unit 
cell or "minimal" coordinate patch defined in ref. |p6|| ). 

Formulating the Hamiltonian on A, with the boundary identifications implied by the 
gauge transformations g(^k), avoids the singularities at the Gribov copies of A = 0. "Bloch 
momenta" associated to the 2vr shift, implemented by the non-trivial homotopy of g(^k), 
label 't Hooft's electric flux quantum numbers |^ "^{Ck = — vr) = exp(vrzefc)\E'(Cfc = vr). 
Note that the phase factor is not arbitrary, but ±1. This is because g'^/^^ is homotopically 
trivial. In other words, the homotopy group of these anti-periodic gauge transformations is 
Z|. Considering a slice of A can obscure some of the topological features. A loop that winds 
around the slice twice is contractible in A as soon as it is allowed to leave the slice. Indeed 
including the lowest modes transverse to this slice will make the Z2 nature of the relevant 
homotopy group evident [jlO|- It should be mentioned that for the torus in the presence 
of fields in the fundamental representation (quarks), only periodic gauge transformations 
are allowed. In that case it is easily seen that the intersection of the fundamental domain 
with the constant abelian gauge fields is given by the domain \Ck\ < 2vr, whose boundary 
coincides with the Gribov horizon. It is interesting to note that points on dQ form an 
explicit example of a continuous degeneracy due to a non-contractible sphere [P^]. 

In weak coupling Liischer showed unambiguously that the wave functionals are 
localized around A = 0, that they are normalizable and that the spectrum is discrete. 
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In this limit the spectrum is insensitive to the boundary identifications (giving rise to a 
degeneracy in the topological quantum numbers). This is manifested by a vanishing electric 
flux energy, defined by the difference in energy of a state with |e| = 1 and the vacuum state 
with e = 0. Although there is no classical potential barrier to achieve this suppression, it 
comes about by a quantum induced barrier, in strength down by two powers of the coupling 



constant. This gives a suppression p8[ with a factor exp{—S/g) instead of the usual factor 
of exp(— STT^/^f^) for instantons Here S = 12.4637 is the action computed from the 
effective potential. At stronger coupling the wave functional spreads out over the vacuum 
valley and the boundary conditions drastically change the spectrum [|TD|. At this point the 
energy of electric flux suddenly switches on. 

The calculation is performed by integrating out the non-zero momentum degrees of 



freedom, for which Bloch degenerate perturbation theory provides a rigorous framework ||3y, 
P], and gives an effective Hamiltonian. Near A = 0, due to the quartic nature of the 
potential energy (F^j)'^ for the zero-momentum modes (the derivatives vanish and the field 
strength is quadratic in the field), there is no separation in time scales between the abelian 
and non-abelian modes. Away from A = one could further reduce the dynamics to one 
along the vacuum valley, but near the origin this would be a singular decomposition (the 
adiabatic approximation breaks down). However, as long as the coupling constant is not 
too large, the wave functional can be reduced to a wave function on the vacuum valley near 
dA where the boundary conditions can be implemented. These boundary conditions are 
formulated in a manner that preserves the invariance under constant gauge transformation 
and the effective Hamiltonian is solved by Rayleigh-Ritz (providing also lower bounds from 
the second moment of the Hamiltonian). The influence of the boundary conditions on the 
low-lying glueball states is felt as soon as the volume is bigger than an inverse scalar 
glueball mass. We summarize below the ingredients that enter the calculations. 

The effective Hamiltonian is expressed in terms of the coordinates c^, where i = {1, 2, 3} 
is the spatial index (cq = 0) and a = {1, 2, 3} is the SU(2)-colour index. These coordinates 
are related to the zero-momentum gauge fields through A°;{x) = cf/L. We note that the 
field strength is given by F^j = —eabdC^C^/L"^ and we introduce the gauge- invariant "radial" 
coordinate = y/J2a ^iCf. The latter will play a crucial role in specifying the boundary 
conditions. For dimensional reasons the effective Hamiltonian is proportional to 1/L. It 
will furthermore depend on L through the renormalized coupling constant {g{L)) at the 
scale n = 1/L. To one-loop order one has (for small L) g{L)^ = 127r^/[— 11 ln(A^f5'L)]. 
One expresses the masses and the size of the finite volume in dimensionless quantities, 
like mass-ratios and the parameter z = mL. In this way, the explicit dependence of g 
on L is irrelevant. This is also the preferred way of comparing results obtained within 
different regularization schemes (i.e. dimensional and lattice regularization) . The effective 
Hamiltonian is now given by 



+7i + 73 1] r^rl + 74 1] rf + 75 ^i^j + 76 H ^ 

i i i>j i i^j i 

+a3 ^-iF^k + «4 E ^'^f + "sdet'c. (10) 

ijk,a ij,a 
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We have organized the terms according to the importance of their contributions, ignor- 
ing terms quartic in the momenta. The first hne gives (when ignoring ai,2) the lowest order 
effective Hamiltonian, whose energy eigenvalues are 0{g'^^^), as can be seen by rescaling 
c with g'^^^. Thus, in a perturbative expansion c = 0{g'^^^). The second line includes the 
vacuum-valley effective potential (i.e. the part that does not vanish on the set of abelian 
configurations). These two lines are sufficient to obtain the mass-ratios to an accuracy of 
better than 5%. The third line gives terms of 0{g^) in the effective potential, that vanish 
along the vacuum-valley. The coefficients (to two-loop order for 7j) are 
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The choice of boundary conditions, associated to each of the irreducible representations 



of the cubic group 0(3, 2L) and to the states that carry electric flux [^, is best described 
by observing that the cubic group is the semidirect product of the group of coordinate 
permutations 5*3 and the group of coordinate reflections Z|. We denote the parity under 
the coordinate reflection c° — c° by Pi = ±1 {i flxed). The electric flux quantum number 
for the same direction will be denoted by gj = ±1. This is related to the more usual 
additive (mod 2) quantum number by qj = exp{inej). Note that for SU(2) electric flux 
is invariant under coordinate reflections. If not all of the electric fluxes are identical, the 
cubic group is broken to S2 x Z|, where 6*2 (~ Z2) corresponds to interchanging the two 
directions with identical electric flux (unequal to the other electric flux). If all the electric 
fluxes are equal, the wave functions are irreducible representations of the cubic group. 
These are the four singlets ^^2)' which are completely (anti-)symmetric with respect to 6*3 
and have each of the parities Pi = ±1. Then there are two doublets E^, also with each of 



(2)- 



Each of these triplet states can 

1(2) 



the parities pi = ±1 and flnally one has four triplets T£ 

be decomposed into eigenstates of the coordinate reflections. Explicitly, for T^^^^ we have 
one state that is (anti-)symmetric under interchanging the two- and three-directions, with 
P2 = P3 = —pi = +1- The other two states are obtained through cyclic permutation of 
the coordinates. Thus, any eigenfunction of the effective Hamiltonian with speciflc electric 
flux quantum numbers qi can be chosen to be an eigenstate of the parity operators pi. The 
boundary conditions of these eigenfunctions ^1^5,^(0) are simply given by 



^9,p(c) 







9 , T 



^f 



PiQi 



-1 



+1 



(12) 



and one easily shows that with these boundary conditions the Hamiltonian is hermitian 
with respect to the innerproduct < \E', \E'' >= (i^c\E'*(c)\E''(c). For negative parity 
states (JliPi = —1) this description is, however, not accurate |]31[ as parity restricted to 
the vacuum valley is equivalent to a Weyl reflection (a remnant of the invariance under 
constant gauge transformations). 
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Figure 4: Mass ratios m' jm as a function of L in units of the inverse scalar glueball 
mass m~^, z = mL. The analytic results are given by the full (continuum) and dashed (4^ 
lattice) curves. Represented are the the square root of the energy of electric flux per unit 
length (bottom), the E tensor mass (middle) and the T2 tensor mass (top). The latter 
two are (almost) degenerate below z = 0.95, which is also where the electric flux energy is 
exponentially suppressed. Where no error bars on the Monte Carlo data are visible they 
are smaller than the size of the data points. 



After correcting for lattice artefacts the (semi-) analytic results agree extremely 
well with the best lattice data |Q (with statistical errors of 2% to 3%) up to a volume of 
about .75 fermi, or about five times the inverse scalar glueball mass. In figure 4 we present 
the comparison for a lattice of spatial size 4^. Monte Carlo data are most accurate for 
this lattice size. For more detailed comparisons see ref. [^. The analytic results below 
z = 0.95 are due to Liischer and Miinster 0, which is where the spectrum is insensitive 
to the identifications at the boundary of A. Apart from the corrections for the lattice 
artefacts, generalization to SU(3) was established by Vohwinkel [Q, with qualitatively 
similar results. In large volumes the rotational symmetry should be restored, as is observed 
from lattice simulations. Most conspicuously the tensor state in finite volumes is split in 
a doublet E, with a mass that is roughly 0.9 times the scalar Ai mass and a triplet T2 
with a mass of roughly 1.7 times the scalar mass. Note that the multiplicity weighted 
average is approximately 1.4 times the scalar mass, agreeing well with what was found at 
large volumes from lattice data ||33|. There has been similar studies using spatial twisted 



boundary conditions (with differing behaviour of the tensor and electric flux states in 



intermediate volumes), see ref. ||35 |. 

At large volumes extra degrees of freedom start to behave non-perturbatively. To 
demonstrate this, the minimal barrier height that separates two vacuum valleys that are 
related by gauge transformations with non-trivial winding number 
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Figure 5: "Artistic impression" of the potential for the three-torus. Shown are two vac- 
uum valleys, related to each other by a gauge transformation hi with winding number 1, 
whose shape in the direction perpendicular to the valley depends on the position along the 
valley. The induced one-loop effective potential, of height 3.210/L, has degenerate minima 
related to each other by the anti-periodic gauge transformations g{^k)- The classical barrier, 
separating the two valleys, has the height 72.605/ Lg"^. 

was found to be 72.605/ Lg"^, using the lattice approximation and carefully taking the 
continuum limit [Q. The situation is sketched in figure 5. One can now easily find for 
which volume the energy of the level that determines the glueball mass (defined by the 
difference with the groundstate energy) starts to be of the order of this barrier height. This 
turns out to be the case for L roughly 5 to 6 times the correlation length set by the scalar 
glueball mass. We expect, as will be shown for the three-sphere, that the boundary of the 
fundamental domain along the path in field space across the barrier (which corresponds to 
the instanton path if we parametrize this path by Euclidean time t), occurs at the saddle 
point (which we call a finite volume sphaleron) in between the two minima. The degrees of 
freedom along this tunnelling path go outside of the space of zero-momentum gauge fields 
and if the energy of a state flows over the barrier, its wave functional will no longer be 
exponentially suppressed below the barrier and will in particular be non-negligible at the 
boundary of the fundamental domain. Boundary identifications in this direction of field 
space now become dynamically important too. The relevant "Bloch momentum" is in this 
case obviously the 6 parameter, as wave functionals pick up a phase factor e*^ under a 
gauge transformation with winding number one. For many of the intricacies in describing 
instantons on a torus we refer to ref. On the three-torus we have therefore achieved 
a self-contained picture of the low-lying glueball spectrum in intermediate volumes from 
first principles with no free parameters, apart from the overall scale. 

4 Gauge fields on the three-sphere 

The reason to consider the three-sphere lies in the fact that the conformal equivalence 
of X IR to IR^ allows one to construct instantons explicitly pSj |T^ , |39[| . This greatly 
simplifies the study of how to formulate 6 dependence in terms of boundary conditions on 
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the fundamental domain, and indeed we will see that for S*^ simple enough results can be 
obtained to address this question The disadvantage of the three-sphere is that in 

large volumes the corrections to the glueball masses are no longer exponential 0. In this 
respect the use of twisted boundary conditions [Bsl offers a viable alternative, but only 



numerical solutions for the relevant instantons are know [pO|] . 

We will summarize the formalism that was developed in |T^. Alternative formula- 
tions, useful for diagonalizing the Faddeev-Popov and fluctuation operators, were given 
in ref. whereas for the explicit formulation of instantons, ref. |^ introduces stere- 
ographic coordinates (demonstrating that simplicity is in the eye of the beholder). We 
embed in IR^ by considering the unit sphere parametrized by a unit vector n^. It is 
particularly useful to introduce the unit quaternions and their conjugates = crj^ by 



a,, 



{id, ir), 



{id, 



-IT] 



(14) 



They satisfy the multiplication rules 



(15) 



where we used the 't Hooft rj symbols p9|, generalized slightly to include a component 
symmetric in ^ and v for a = 0. We can use 77 and f] to define orthonormal framings of 
S^, which were motivated by the particularly simple form of the instanton vector potentials 
in these framings. The framing for S"^ is obtained from the framing of by restricting 
in the following equation the four-index a to a three-index a (for a = one obtains the 
normal on 5''^): 

e° = t^u^v, e° = fi'^^n^. (16) 

Note that e and e have opposite orientations. Each framing defines a differential operator 
and associated (mutually commuting) angular momentum operators Li and L2: 



d 



(17) 



It is easily seen that L\ = L\, which has eigenvalues /(/ + 1), with / 

The (anti-) instantons |Q in these framings, obtained from those on IR^ by interpreting 
the radius in IR^ as the exponential of the time t in the geometry S'^ x IR, become 



Oil 

)4 



An 



e ■ a 



2{l + e-ny 



A 



a A e— {u + e ■ n)a 
2{l + e-n) ' 



where 



2s^ 



u 



2sb„ 



1 + 62 + s2 ^ 



Ae*. 



(19) 



Here e and A are defined with respect to the framing for instantons and with respect 
to the framing eJJ for anti-instantons. The instanton describes tunnelling from A = at 
t = —00 to Aa = —cTa Sit t = oo, over a potential barrier at t = that is lowest when = 0. 
This configuration corresponds to a sphaleron |Q, i.e. the vector potential Aa = —\(Ja is 
a saddle point of the energy functional with one unstable mode, corresponding to the 
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direction (u) of tunnelling. At t = oo, Aa = —era has zero energy and is a gauge copy of 
= by a gauge transformation g = n ■ a with winding number one. 
We will be concentrating our attention to the modes that are degenerate in energy to 
lowest order with the modes that describe tunnelling through the sphaleron and "anti- 
sphaleron" . The latter is a gauge copy by a gauge transformation g = n ■ a with winding 
number —1 of the sphaleron. The two dimensional space containing the tunnelling paths 
through these sphalerons is consequently parametrized by u and v through 

A,{u,v) = {-ue;,-ve;)^. (20) 

The gauge transformation with winding number —1 is easily seen to map {u,v) = {w,0) 
into {u, v) = {0,2 — w). The 18 dimensional space is defined by 

A,{c, d) = (c^; + d^ej,) ^ = Aic d)el. (21) 

The c and d modes are mutually orthogonal and satisfy the Coulomb gauge condition: 

dMc,d) = Q. (22) 

This space contains the {u,v) plane through = —u5^ and d1 = —vS^. The significance 
of this 18 dimensional space is that the energy functional |T^ 



V(c, d)^- J^^ \ tr (4) = V(c) + V{d) + ^ {(cn^(4)^ - [c'^d'^f] , (23) 

V(c) = 27r^ [2{dtf + Qdetc+\[{cXf - [c^c'^f]] , (24) 

is degenerate to second order in c and d. Indeed, the quadratic fluctuation operator ^A in 
the Coulomb gauge, defined by 

-/^^itr(4) = /^^tr(A,MA) + (^(^'), 

M., = 2L% + 2(L, + Sy,,, S^ = -zea,j, (25) 



has A{c, d) as its eigenspace for the (lowest) eigenvalue 4. These modes are consequently the 
equivalent of the zero-momentum modes on the torus, with the difference that their zero- 
point frequency does not vanish. An effective Hamiltonian for the c and d modes is here 
derived from the one-loop effective action and errors due to an adiabatic approximation are 
not necessarily suppressed by powers of the coupling constant. Nevertheless, one expects to 
achieve an approximate understanding of the non-perturbative dynamics in this way 



FP\[A) in eq. (^ is defined as an hermitian operator acting on the vector space C of 
functions g over with values in the space of the quaternions IH = {q^(y^\q^ G IR}. The 
gauge group Q is contained in C by restricting to the unit quaternions: Q = {g & C\g = 
g^iO'fj,,gii G 'Si^g^g^ = !}• For arbitrary gauge groups C is defined as the algebra generated 
by the identity and the (anti-hermitian) generators of the algebra. When minimizing the 
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u 

Figure 6: Location of the classical vacua (large dots), sphalerons (smaller dots), the 
Gribov horizon (fat sections), the boundary of A (dashed curves) and part of the boundary 
of the fundamental domain (full curves). Also indicated are the lines of equal potential in 
units of 2" times the sphaleron energy. . 



same functional over the larger space L one obviously should find a smaller space A C A. 
Since £ is a linear space A can also be specified by the condition that FP\ (A) be positive, 

A = {AGr|((?,FPi(A)(7)>0, V(7G/:}. (26) 

As for the Gribov horizon, the boundary of A is therefore determined by the location where 
the lowest eigenvalue vanishes. For the (c, ci) space it can be shown ||25| that the boundary 
dk will touch the Gribov horizon 9^2. This establishes the existence of singular points on 
the boundary of the fundamental domain due to the inclusion A C A C fi. By showing 
that the fourth order term in eq. (^ is positive (see app. A of ref. p5[) this is seen to 
correspond to the situation as sketched in figure 1. 

One can make convenient use of the SU(2)^ symmetry generated by Li, L2 and T to 
calculate explicitly the spectrum of FPt^A). One has 

FP,(A(c, d)) = ALj - -^ctT^V, - -^dtT^L\, (27) 

which commutes with L\ = but for arbitrary {c,d) there are in general no other 
commuting operators (except for a charge conjugation symmetry when t = \). Restricting 
to the (m, v) plane one easily finds that 

FPt{A{u, v)) = AL\ + juU ■ ft + jvU ■ ft, (28) 

which also commutes with the total angular momentum = Li + L2 + Tt and is easily 
diagonahzed. Figure 6 summarizes the results for this (m, v) plane and also shows the 
equal-potential lines as well as exhibiting the multiple vacua and the sphalerons. As it 
is easily seen that the two sphalerons are gauge copies (by a unit winding number gauge 
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X3 




X3 




Figure 7: The fundamental domain (left) for constant gauge fields on S^, with respect 
to the "instanton" framing e^, in the diagonal representation Aa = Xa^a (no sum over a). 
By the dots on the faces we indicate the sphalerons, whereas the dashed lines represent 
the symmetry axes of the tetrahedron. To the right we display the Gribov horizon, which 
encloses the fundamental domain, coinciding with it at the singular boundary points along 
the edges of the tetrahedron. 



transformation) with equal norm, they lie on (9A, which can be extend by perturbing around 
these sphalerons [Q. 

To obtain the result for general (c, d) one can use the invariance under rotations gener- 
ated by Li and L2 and under constant gauge transformations generated by Tt, to bring c 
and (i to a standard form, or express det [F Pt{A{c, d)) 1 ^ , which determines the locations 

of dVt and (9A, in terms of invariants. We define the matrices X and Y by = (cc*)JJ and 
= {dd^)l, which allows us to find 

det (f Pi {A{c,d))\i^i) = [81 - 18 Tr (X + F) + 24(det c + det d) 



tt{x -Y)y + 2Tt{{x -Yy)Y. (29) 



The two-fold multiplicity is due to charge conjugation symmetry. The expression for t = 1, 
that determines the location of the Gribov horizon in the (c, d) space, is given in app. B 



of ref. If we restrict to d = the result simplifies considerably. In that case one 

can bring c to a diagonal form = Xi6f. The rotational and gauge symmetry reduce to 
permutations of the Xi and simultaneous changes of the sign of two of the Xi. One easily 
finds the invariant expression ( Tr (X) = J2i ^1 and det c = Hi ^j) 

det (fPi(A(c,0))|;^i) = (2detc-3Tr(X) + 27)\ (30) 

In figure 7 we present the results for A and Vt. In this particular case, where c? = 0, 
A coincides with A, a consequence of the convexity and the fact that both the sphalerons 
(indicated by the dots) and the edges of the tetrahedron lie on 9A, the latter also lying 
on dVL. It is essential that the sphalerons do not lie on the Gribov horizon and that 
the potential energy near dVl is relatively high. This is why we can take the boundary 
identifications near the sphalerons into account without having to worry about singular 
boundary points, as long as the energies of the low-lying states will be not much higher 
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than the energy of the sphaleron. It allows one to study the glueball spectrum as a function 
of the CP violating angle 9, but more importantly it incorporates for ^ = the noticeable 
influence of the barrier crossings, i.e. of the instantons. For details see . 



5 Discussion 

We have analysed in detail the boundary of the fundamental domain for SU(2) gauge 
theories on the three-torus and three-sphere. It is important to note that it is necessary to 
divide A by the set of all gauge transformations, including those that are homotopically 
non-trivial, to get the physical configuration space. All the non-trivial topology is then 
retrieved by the identifications of points on the boundary of the fundamental domain. 

As we already mentioned in the introduction, the knowledge of the boundary identi- 
fications is important in the case that the wave functionals spread out in configuration 
space to such an extent that they become sensitive to these identifications. This happens 
at large volumes, whereas at very small volumes the wave functional is localized around 
y4 = and one need not worry about these non-perturbative effects. That these effects 
can be dramatic, even at relatively small volumes (above a tenth of a fermi across), was 
demonstrated for the case of the torus. However, for that case the structure of the funda- 
mental domain (restricted to the abelian zero-energy modes) is a hypercube and deviates 
considerably from the fundamental domain of the three-sphere. Results for the spectrum 



in the latter case recently became available |^ and indicate that the tensor to scalar glue- 



ball mass ratio is compatible in volumes around one fermi. For the three-sphere the tensor 
glueball is of course not split into a doublet and triplet representation. 

It should be noted that the shape of A is independent of L if the gauge field is expressed 
in units of 1/L. Suppose that the coupling constant will grow without bound. This would 
make the potential irrelevant and makes the wave functional spread out over the whole of 
field space (which could be seen as a strong coupling expansion). If the kinetic term would 
have been trivial the wave functionals would be "plane waves" on a space with complicated 
boundary conditions. In that case it seems unavoidable that the infinite volume limit would 
depend on the geometry (like or S'^) that is scaled-up to infinity. With the non-triviality 
of the kinetic term this conclusion cannot be readily made and our present understanding 
only allows comparison in volumes around one cubic fermi. However, one way to avoid this 
undesirable dependence on the geometry is that the vacuum is unstable against domain 
formation. As periodic subdivisions are space filling on a torus, this seems to be the 
preferred geometry to study domain formation. In a naive way it will give the correct 
string tension (flux conservation tells us to "string" the domains that carry electric flux) 
and tensor to scalar mass ratio (averaging over the orientations of the domains is expected 
to lead to a multiplicity weighted average of the T2 and E masses). Furthermore, the 
natural dislocations of such a domain picture are gauge dislocations. The point-like gauge 
dislocations in four dimensions are instantons and in three dimensions they are monopoles. 
Their density is expected to be given roughly as one per domain (with a volume of around 
0.5 cubic fermi). Also the coupling constant will stop running at the scale of the domain 



size. We have discussed this elsewhere and refer the reader to refs. [2(:, it] for further 



details, as the ideas in this direction remain speculative. In the context of twisted boundary 



conditions, related ideas were recently developed in ref. [47 . 
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